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Abstract. We have re-analysed all available high-resolution ultraviolet IUE spectra of the high-mass X-ray binary 
HD153919/4U1700-37. The radial velocity semi-amplitude of 20.6 ± 1.0 km s _1 and orbital eccentricity of 0.22 ± 0.04 
agree very well with the values obtained earlier from optical spectra. They disagree with earlier conclusions for the same data 
reduced by Heap & Corcoran 1 1992 1 and by Stickland & Lloyd 1 1993 1. 
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1. Introduction 

The bright Of star HD153919 is the optical counterpart of the 
3.4 day eclipsing high-mass X-ray binary 4U 1700-37, discov- 
ered by Jones et al. dl973> . About 2.5 million years ago, the 
binary system escaped from the OB association Sco OB2 after 
the supernova explosion that formed the compact companion 
(Ankay et al. 2001). The X-ray source powered by the stellar 
wind of the Of star shows no regular pulsations and the highly 
variable X-ray flux is less than expected from standard wind 
accretion (e.g. Rubin et al. 1996 Kaper 1998). The decrease 
in orbital period determined from mid-eclipse observations and 
first mentioned by Haberl et al. ( 1989 ) was confirmed by Rubin 
et al. 



1199 6 1 and is probably due to wind-driven angular mo- 
mentum loss from the Of star. Because no direct measurements 
can be made of the orbital motion of the compact object, sys- 
tem parameters depend heavily on radial velocity studies of the 
optical companion. The most recent analysis of the physical 
parameters cannot resolve the nature of the compact object in 
this system (Clark et al. 2002 1. Its likely mass is considerably 
lower than those found for black holes, but it is significantly in 
excess to fit the high density nuclear equation of state for neu- 
tron stars. A point to remember, however, is that also the Vela 
X-l pulsar has a high mass, comparable to that of 4U 1700-37 
(van Kerkwijk et al. 119951 . 

The most thorough study of radial velocities of lines in 
the visible spectral region has been made by Hammerschlag- 
Hensberge (1978 1. It is clear from the observations that the 



strongly variable stellar wind in this system complicates a 
straightforward determination of the system parameters. The 
results for the orbital analysis of all lines together in the visible 
blue spectral region (3700 — 4900 A) yield a semi-amplitude 
K = 19±lkms _1 and an orbital eccentricity e = 0.16±0.08. 



Send offprint requests to: G. Hammerschlag-Hensberge; e-mail: 
godeliev @ science.uva.nl 

* Based on observations obtained with the International Ultraviolet 
Explorer 



Heap & Corcoran ( 1992 1 have examined the system param- 
eters using the available ultraviolet spectra of HD153919. They 
find K = 18 ± 3 km s _1 and a negligible eccentricity. For 
their investigation they cross-correlated a few selected narrow 
regions of the ultraviolet spectra with the same regions of one 
of the spectra. To bring their spectra to the same frame of ref- 
erence they used the heliocentric correction. However, instru- 
mental and processing shifts during the reduction of spectra 
obtained with the International Ultraviolet Explorer (IUE) re- 
quire a more sophisticated approach (Stickland 1992 1. 

Stickland & Lloyd ( 1993 1 re-examined the same data, now 
cross-correlating both stellar and interstellar spectra with refer- 
ence to one of the UV spectra. They used the complete ultravi- 
olet spectral range excluding P-Cygni lines whose profiles are 
affected by mass loss. Amazingly, they found a much smaller 
semi-amplitude K of 10 km s _1 and an eccentricity of 0.17. 

In a similar study of the high-mass X-ray binary 
HD77581/Vela X-l Barziv et al. (120011 found a disagree- 
ment in the semi-amplitude of the radial velocity curve with 
the results of Stickland et al. ( 1997 ) for this system. This 
strengthened our decision to re-analyse the ultraviolet spectra 
of HD153919 again. In Sect. 2 we describe the observations 
and the data reduction procedures. In Sect. 3 we present the re- 
sulting radial velocity curve and in Sect. 4 we compare our re- 
sults with the ones found by Heap & Corcoran and by Stickland 
& Lloyd. In Sect. 5 we summarize our conclusions and discuss 
the evolutionary consequences. 
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Table 1. Journal of IUE SWP high-resolution spectra of HD153919. Image numbers include an S or an L, refering to the small and 
large aperture, respectively. Phases are calculated from the ephemeris T = JD2448900.873 ± 0.015; P = 3.411581 ±0.000027 
given by Rubin et al. ( 1996 1 including the orbital period change P/P — —(3.3 ± 0.6) x 10~ 6 yr _1 . The velocities are corrected 
for a shift of interstellar lines and are relative to that of SWP 1476 which was arbitrarily set to zero. 
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2. Observations and data reduction 



A large number of ultraviolet spectra of HD153919 have been 
obtained with the International Ultraviolet Explorer. The spec- 
tra obtained with the Short Wavelength Prime Camera in high- 
resolution mode are listed in Tabled A separate analysis of the 
same IUE spectra was performed by Stickland & Lloyd ( 1993 1. 
These authors cross-correlated the set of spectra with one par- 
ticular spectrum (SWP25600). Because their results were in 
large discrepancy with results from optical data we decided to 
re-analyse the same IUE data. We cross-correlated each spec- 
trum with all other spectra of the data set to reduce errors. 
We did not use a template spectrum of another star with iden- 
tical spectral type as independent reference spectrum- this is 
not realistic for such hot mass-losing stars which each have 
unique spectra- nor did we use an average spectrum to avoid 
auto-correlation effects leading to a systematic decrease of the 
amplitude of the radial velocity curve. 



First the spectra - initially mapped onto a uniform wave- 
length grid of 0.05 A - were transformed to a logarithmic wave- 
length scale such that the Doppler shift becomes a linear dis- 
placement along the spectrum. The relevant part of the cross- 
correlation function is the area close to the top; it is this up- 
per part which is fitted with an analytic function (a Gaussian 
plus a linear function) whose maximum determines the veloc- 
ity shift of the spectrum. Because the cross-correlation pro- 
files of IUE spectra show small excess peaks (see Evans ( 1988 1 
and Fig. 4 in Van Kerkwijk et al. ( 1995)) due to the so-called 
fixed-pattern noise in the detector, we corrected for those by 
including an extra Gaussian peak in the fit. For details on the 
cross-correlation method and instructive figures we refer to Van 
Kerkwijk ( fT993t . 

The spectral regions which were used for the cross- 
correlation are listed in Table |2] The regions were chosen in 
a way to exclude parts of the spectrum containing variable line 
profiles such as the well-known P-Cygni profiles (see for in- 
stance Hammerschlag-Hensberge et al. ( 1990)); also the spec- 
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Table 2. Wavelength regions used for cross-correlation. 
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tral region including variable emission lines 1540 — 1740 A 
(Kaper et al. 1990 1 was excluded. Table|2]lists the wavelength 
regions including interstellar lines which were used to place 
the individual spectra onto the same frame of reference, be- 
cause IUE spectra have a variety of instrumental and process- 
ing shifts, as mentioned before. 

3. Results 

In Table [Qwe list the radial velocities derived from the cross- 
correlations of all spectra with each other. In Fig. Q we 
show the radial velocities we obtained from the IUE spec- 
tra as a function of orbital phase. Because the spectra were 
cross-correlated with each other, no 7-velocity has been deter- 
mined for the system and the velocities in the figure have been 
shifted to zero system velocity. The 7-velocity has been deter- 
mined elsewhere (e.g. Gies & Bolton 1986 Humphreys 1978 1: 
—60 km s _1 . It was not included in our figure because the value 
- which is rather uncertain - did not follow from our data. At 
first glance, it is already clear that the data show significant de- 
viations from a smooth Keplerian curve. We believe these are 
intrinsic to the star, since earlier optical observations have al- 
ready shown that the star has intrinsic variations (e.g., Kaper et 
al.[H>94i. As we will find that the deviations are substantially 
larger than the measurement errors, we give each measurement 
equal weight in our fits below. 

We first fit a circular orbit, shown by the long-dashed line 
in Fig.Q an d find K = 18-7 ± 1.0 km s _1 . The root-mean- 
square (rms) residual is 6.36 km s _1 . This is a factor two larger 
than the mean error in the measured velocities (using the formal 
errors, one finds a reduced \ 2 — 3.6). 

The residuals relative to the best-fit circular orbit are shown 
in the lower panel of Fig. Q These seem to show a double- 
wave pattern, suggestive of an eccentric orbit. Fitting an ec- 
centric orbit, we find eccentricity e = 0.22 ± 0.04, perias- 
tron angle uj = 49° ± 11°, and the radial-velocity amplitude is 
K = 20.6 ± 1.0 km s -1 . The rms residuals are slightly lower, 
at 5.85 km s _1 . Doing an F-test, this is not a significant im- 
provement. The F-test, however, ignores the correlations in the 
residuals. In order to obtain a more robust estimate of the sig- 
nificance, we ran Monte-Carlo simulations, in which we made 
artificial data sets in which the velocity at each measurement 
phase was given by the sum of the velocity expected from the 
best-fit circular orbit and a random deviation. For the latter, 
we used two components, one to mimic the measurement error 
(i.e., a random number drawn from a normal distribution with 
standard deviation equal to the measurement error), and one to 
mimic the intrinsic deviations. For the intrinsic deviations, we 
again assumed a normal distribution, with a standard deviation 



of 6 km s _1 , so that our artificial data sets have the same rms 
residuals as our data. Our simulations confirmed our error es- 
timate on the eccentricity, and showed that the probability of 
obtaining an eccentricity of 0.22 by chance due to the devia- 
tions was negligible. 

We should note that what is proven above is that for a circu- 
lar orbit there are systematic variations with orbital phase over 
and above the significant random excursions. An eccentric orbit 
is a simple explanation for these systematic variations, but we 
cannot exclude other systematic effects. A cause for caution is 
that Barziv et al. (2001 1 found from their optical data an eccen- 
tricity which was larger than the value from the X-ray pulsar in 
VelaX-1 (see their Fig. 18). Furthermore, model calculations of 
Zuiderwijk et al. (I1977> showed that the radial-velocity curve 
for a tidally distorted star in a circular orbit might show an ap- 
parent eccentricity. The shape of the expected curve, however, 
is different from that observed. Given this, we cannot dismiss 
the possibility of a circular orbit out of hand, although we be- 
lieve that the simplest and most likely solution is that it is in 
fact eccentric. 

It is obvious that for some spectra the individual radial- 
velocity measurements deviate significantly from the best-fit 
Keplerian orbit. This was already known from the optical spec- 
tra. The Of supergiant has a strong stellar wind which can cause 
irregular variations in the radial velocities of some lines. Also 
the presence of a photo-ionization wake (Kaper et al. 0994 1 
trailing the X-ray source may cause irregular small orbit-to- 
orbit variations in the radial velocities. Nevertheless we de- 
voted extra attention to some of the spectra. Particularly strik- 
ing are the three velocities near phase 0.25 that are much 
lower than the other spectra at the same phase interval. These 
three are for spectra SWP1982, 1983 and 1986, which were 
taken in sequence. This would suggest an instrumental prob- 
lem, were it not that SWP1987, which was taken immedi- 
ately after SWP1986 at the same night as the mentioned three 
spectra, has a much higher velocity comparable with that of 
SWP 4751, 4752 and 4753. We compared the line profiles of 
all those spectra but could not find significant differences, al- 
though SWP 1982 and 1983 are very noisy. Note that the mod- 
est spectral resolution (R ~ 10, 000) and low S/N (~ 20) does 
not allow for a detailed study of individual photospheric lines. 

For HD 77581, the companion of Vela X-l, it was found 
that sequences of radial velocities showed deviations that 
were correlated over periods of about one day (Van Kerkwijk 
et al. I1995> . To see whether such deviations are present in 
HD 153919 as well, in Fig.^lines connect residual velocities 
for spectra that were taken as part of a sequence. Unfortunately, 
the result is inconclusive. Some groups seem to show correlated 
deviations, but for others the deviations seem random. It would 
be certainly worthwhile to monitor HD 153919 during some 
nights continuously to have a clue to the origin of those ex- 
cursions and stellar wind variations that certainly influence and 
hamper an accurate determination of the orbital parameters. 

Finally, we should mention that almost half of the set of 
spectra was taken in the small aperture of the SWP camera, the 
others in the large aperture (see TableJQ. The spectra taken in 
the small aperture are the ones of the first years of IUE and 
have generally lower S/N. We determined orbital solutions for 
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Fig. 1. Radial-velocity curve derived from the cross-correlated IUE spectra of HD153919. The radial velocities in the top panel 
are plotted with error bars indicating the la uncertainties. Overdrawn is the Keplerian curve that best fits the data (short dashes) 
as well as the solution for a circular orbit (long dashes). The velocities have been shifted to zero system velocity based on 
the best fitting circular orbit (our analysis did not allow us to derive the systemic velocity; see text). In the lower panel, the 
residuals relative to the best-fit circular orbit are shown; the short-dashed curve shows where the residuals should lie for the 
best-fit eccentric orbit. For clarity, the error bars have been omitted. Points connected by lines were taken as a sequence. 



the small and large-aperture sets separately, but did not find 
significant differences. 

4. Comparison with earlier results 

From the present results we can conclude that the derived radial 
velocities are in good agreement with the results obtained much 
earlier from optical data (see Hammerschlag-Hensberge 1978 1. 

Our results do not differ much from those by Heap & 
Corcoran ( 1992 1 for the same data. The advantage of our anal- 
ysis is that we use a larger part of the spectrum and correlate all 
spectra with each other, which, as a consequence, results in a 
smaller error on K. Remarkably, Heap & Corcoran ( 1992 1 find 
a negligible eccentricity (e ~ 0.003) with a very small error of 
the order of 0.003. At first glance, the radial velocities in their 
Fig. 2 seem to fit equally well to an eccentric orbit. Also, from 
a statistical point of view, their errors on e and u> are much too 



small. For small eccentricities the deviation from a circular or- 
bit is a sinusoid with half the binary period and it is easy to see 
that e ~ A(P/2) /A(P) = A(P/2) /K, where A is the ampli- 
tude of the sinusoid and K the radial velocity amplitude. The 
error on A(P/2) will be comparable to the one on K and as 
a consequence a e ~ <tk/K . This implies er e ~ 0.17 for the 
error on the eccentricity, instead of <r e = 0.003 found by Heap 
& Corcoran. 

On the other hand, Stickland & Lloyd J1993I find a stellar 
orbit with an eccentricity close to ours, but their results show 
a much smaller semi-amplitude K of 10 km s _1 . Contact by 
one of the authors with Dave Stickland for similar discrepant 
velocity values in the X-ray binary Vela X-l has proven that the 
main cause for the discrepancy in the results is that their cross- 
correlation procedure sets to zero the regions around interstel- 
lar lines and wind lines. When this is combined with slight in- 
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adequacies in normalization, artificial lines are introduced at 
fixed wavelengths in both spectra to be cross-correlated. As 
a consequence the radial velocities are biased systematically 
towards zero and the radial-velocity amplitude will be under- 
estimated. A discrepancy of about 4 km s _1 had been found 
between the radial-velocity amplitude for the ultraviolet spec- 
tra of Vela X-l analyzed by both groups (see Barziv et al. 2001 
and Stickland et al. H997i . The reason for the much larger dis- 
crepancy of 10 km s _1 for HD153919 might be the weaker and 
broader photospheric lines in this star combined with the bet- 
ter signal-to-noise ratio which give more striking errors when 
interstellar and wind lines are set to zero. 

5. Conclusions and evolutionary consequences 

The two main conclusions of this paper are that (i) the radial- 
velocity amplitude K of HD153919 determined from ultravi- 
olet spectra is consistent with that derived from optical spec- 
tra; and (ii) the orbit most likely has a significant eccentric- 
ity. Note that the obtained value for K depends on whether 
the orbit is eccentric (K = 20.6 km s~\ e = 0.22) or not 
(K = 18.7 km s~ l ). 

Clark et al. ( 2002 1 performed a thorough study of the sys- 
tem parameters of HD153919/4U1700-37 and determined the 
most likely masses of both components using Monte Carlo sim- 
ulations. These suggest a mass for HD153919 and 4U1700-37 
of M* = 58 ± 11 M Q and M x = 2.4 ± 0.3 M Q , respectively. 
Clark et al. used the results described in this paper in their anal- 
ysis, but prefered the orbital solution with zero eccentricity. As 
the mass function depends on the eccentricity: 

/ = 1.038 x 1Q- 7 K 3 P(1 - e 2 ) 3/2 (1) 

and 

, 9 (M x sini) 3 
(M* + Mx) = V X 1 , (2) 

they noticed that the masses of both components become sig- 
nificantly higher when taking e = 0.22. They argue that such 
high values for M* are inconsistent with log g derived from the 
analysis of the photospheric spectrum. A non-zero eccentricity 
also results in a higher number of solutions that are rejected 
due to the inclination constraints. However, Clark et al. took 
the value of the radial-velocity amplitude belonging to the ec- 
centric orbit (K = 20.6 km s _1 , so that / = 0.0031 M Q ), 
rather than K = 18.7 km s -1 for the circular orbit (/ = 
0.0023 M©); for the latter value the masses would have be- 
come even higher. If one takes K = 20.6 km s _1 and e = 0.22 
then / = 0.0029 Mq. Therefore, we conclude that the masses 
of both components are about 4 % higher than proposed by 
Clark et al. ( 2002 1 if one uses the correct value for / in case of 
the eccentric orbital solution. 

The short orbital period of the system, combined with its 
eccentricity, will induce strong tidal interactions between the 
two components. Quaintrell et al. (2003) present evidence for 
tidally induced non-radial oscillations in HD77581, the B- 
supergiant companion to Vela X-l, like 4U 1700-37 an eccen- 
tric (e = 0.09), wind-fed high-mass X-ray binary with an 
orbital period of 8.9 days. These non-radial oscillations may 



well be the cause of the radial-velocity excursions detected in 
Vela X-l. The radial-velocity excursions in HD153919 could 
have a similar origin. 

The eccentricity of the orbit of both Vela X-l and 4U1700- 
37 suggests that these systems have just entered the high-mass 
X-ray binary phase and are on their way to become Roche- 
lobe overflow systems. The latter systems have circular or- 
bits, possibly the result of tidal interaction. Both HD77581 
and HD153919 do not rotate synchronously with their orbit: 
-Porb/-frot is 0.67 (Zuiderwijk 1995) and 0.46 (taking v sin i = 
150 km s~ 1 , i = 90°, i?* = 21.9 R©), respectively, which may 
indicate that the supergiants are still expanding and that tidal 
forces have not (yet?) managed to synchronize the system and 
circularize the orbit. 

Conti ( 1978 1 and Rappaport & Joss ( 1983) suggested that 
the OB supergiants in HMXBs are too luminous for their mass. 
Kaper (2001) showed that the discrepancy between the mea- 
sured and "spectroscopic" mass is most severe for the Roche- 
lobe overflow systems, but is only modest for wind-fed sys- 
tems such as Vela X-l. The proposed high mass of HD153919 
is consistent with its high luminosity, again suggesting that 
HD153919 has not yet reached the phase of Roche-lobe over- 
flow. 
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